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^ . Abstract 
. 

We report on absolute proper-motion measurements of an H2O maser source in the NGC 281 West 
molecular cloud, which is located ~320 pc above the Galactic plane and is associated with an Hi loop 
^ ■ extending from the Galactic plane. We have conducted multi-epoch phase-referencing observations of the 

ly^ ' maser source with VERA (VLBI Exploration of Radio Astrometry) over a monitoring period of 6 months 

since May 2006. We find that the H2O maser features in NGC 281 West are systematically moving toward 
CO ' the southwest and further away from the Galactic plane with a vertical velocity of ^^20— 30 km s^^ at its 

, estimated distance of 2.2—3.5 kpc. Our new results provide the most direct evidence that the gas in the 

' NGC 281 region on the H I loop was blown out from the Galactic plane, most likely in a superbubble driven 

' by multiple or sequential supernova explosions in the Galactic plane. 

Key words: Galaxy: kinematics and dynamics — ISM: bubbles — ISM: H II regions — ISM: individual 
(NGC 281) — masers (H2O) 



■ 1. Introduction considered to have been formed by multiple or sequential 

JH \ supernovae and stellar winds from OB associations (e.g.. 

Over the past few decades, numerous efforts have been Tomisaka et al. 1981; Tomisaka &: Ikeuchi 1986; Mac Low 

made both observationally and theoretically to under- et al. 1989; Tomisaka 1998). An excellent review of the 

stand the large-scale superstructures such as superbub- field is given by Tenorio-Tagle and Bodenheimer (1988). 

bles/supershells in the interstellar medium in the Galaxy Subsequent studies of the superstructures have led to an 

and in nearby galaxies. These large-scale superstructures extended study of the disk-halo interaction in the Galaxy 

are characterized by their circular or arc-like appearances and in other galaxies (Mac Low et al. 1989; Norman & 

with radii ranging from tens of parsecs up to 1000 pc. Ikeuchi 1989) and, as represented by the concept of galac- 

Many of these Galactic shells were found by Heiles (1979) tic 'chimneys', the blowout mechanism of superbubbles 

based on the survey of the Hi 21-cm line by Weaver and has drawn attention as a form of efficient energy transport 

Williams (1973). The energy required for such a large- from a galaxy's disk to its halo (e.g., Norman & Ikeuchi 

scale superstructure with a radius greater than 100 pc 1989; Normandeau et al. 1996). 

ranges up to 10^^ erg and is hundreds of times larger than For studying the origin, energetics and the timescale of 
that available from a single supernova explosion (Heiles superbubbles/supershells, it is of great importance to in- 
1979). For many years, the long-term evolution of super- vestigate the kinematics and dynamics of atomic or molcc- 
nova remnants (SNRs) has been studied (e.g.. Cox 1972; ular clouds and Hll regions associated with these super- 
Chevalier 1974; Tomisaka & Ikeuchi 1986; Hanayama & structures. Whereas theoretical simulation models of su- 
Tomisaka 2006), and those superbubbles/supershells are pernova evolution have made considerable progress and 
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become increasingly more complex, there had been hardly 
any means to measure the proper motions of superstruc- 
ture blowouts obscrvationally, due to insufRcicnt astro- 
metric accuracy of available instruments. However, recent 
developments in VLBI technology, especially the phase- 
referencing VLBI technique, provide a powerful tool with 
its highest possible precision for absolute astrometry to 
overcome these observational difficulties and offer a break- 
through in our understanding of the superstructure mech- 
anism in the Galaxy. 

The NGC 281 (Sh 184; at ajooo =00'^52'", S2000 = 
-H56°34' or I =123°. 07, b = -6.°31) region is a good ob- 
servational target for studying the blowout dynamics of 
superbubbles. Based on Hi data from Hartmann and 
Burton (1997), Megeath et al. (2002) and Megeath et al. 
(2003) identified the NGC 281 molecular cloud complex, 
situated z ~320 pc above the midplane of the Perseus arm 
of the Galaxy, on a large-scale H I loop from the Galactic 
plane. They suggested that these molecular clouds at such 
a remarkable vertical height from the Galactic plane have 
formed in the gas swept up and compressed in a blowout 
triggered by multiple supernova explosions. 

The NGC 281 nebula is an Hll region surrounded by a 
giant molecular cloud complex which was entirely mapped 
in the CO emission line by Lee and Jung (2003). The de- 
termination of the photometric/kinematic distance to the 
NGC 281 region has been attempted for many years, but 
the obtained values vary widely from 2 kpc to 3.5 kpc 
(e.g., Sharpless 1954; Georgelin & Georgelin 1976; Roger 
& Pedlar 1981; Henning et al. 1994; Guetter & Turner 
1997; Lee & Jung 2003). The source of ionization for the 
NGC 281 nebula is a star cluster whose brightest mem- 
ber is a compact, Trapezium-like OB multiple star system 
caUcd ADS 719 or HD 5005 (Sharpless 1954; Elmegreen & 
Lada 1978; Henning et al. 1994; Guetter & Turner 1997). 
A southwestward proper motion of HD 5005 was observed 
with Hipparcos (Perryman et al. 1997), but with insuffi- 
cient precision to discuss the motion of the system. 

Because of the clear separation between the H II region 
and the surrounding molecular clouds, along with its high 
Galactic latitude with negligible contamination by inter- 
stellar matter, the NGC 281 region has been a good labo- 
ratory for studying star formation at multiple wavelengths 
for more than a decade. Elmegreen and Lada (1978) dis- 
covered an H2O maser source at the ^^CO emission peak 
in the NGC 281 West molecular cloud (located ~5' to 
the southwest of HD 5005) near the molccular/HlI in- 
terface and proposed that the expansion of the Hll re- 
gion is triggering star formation in the adjoining molec- 
ular cloud. At this location of the ^^CO emission peak 
in the NGC 281 West cloud, there is a far-infrared source 
IRAS 00494-1-5617, and an embedded cluster of low mass 
stars has also been detected in the near-infrared obser- 
vations (Carpenter et al. 1993; Hodapp 1994; Megeath 
& Wilson 1997). Tofani ct al. (1995) observed the H2O 
maser source in NGC 281 West cloud with the VLA and 
found that the H2O maser emission is resolved into three 
distinct spatial components CI, C2 and C3 with maxi- 
mum separation of ~8". 



The H2O maser source in the NGC 281 West cloud pro- 
vides an excellent target for high-resolution VLBI obser- 
vations. The average of the maser motions is expected to 
trace the systemic motion of the molecular cloud on the 
H I loop, which then yields a direct measure of the super- 
bubble expansion. In order to detect the systemic motion 
of the NGC 281 West cloud on the Hi loop, we therefore 
conducted high-resolution observations of the H2O maser 
source in the cloud with VERA, a Japanese VLBI ar- 
ray dedicated to phase- referencing VLBI astrometry (e.g., 
Honma et al. 2000; Kobayashi et al. 2003; Honma ct al. 
2005; Imai ct al. 2006). In this paper, wc report on our 
successful detection of the absolute proper motion of the 
H2O maser source, which reveals the motion of this re- 
gion away from the Galactic plane, most likely having 
originated in the Galactic plane and been blown out by a 
superbubblc blowout. 

2. VERA Observations and Data Reduction 

The VERA observations were conducted at six epochs, 
spaced approximately at monthly intervals, on 2006 May 
14, July 21, August 3, September 5, October 25, and 
November 18 (day of year 134, 202, 215, 248, 298, 
and 322, respectively). In each epoch, the observation 
was carried out for 7—9 hours in the dual-beam mode 
for phase referencing (e.g., Honma et al. 2003). The 
calibration data of real-time phase difference between 
the two beams were taken during the observations, us- 
ing artificial noise sources injected into two beams at 
each station (Kawaguchi et al. 2000). An H2O maser 
source in the NGC 281 West cloud was observed si- 
multaneously with one of the two alternately-scanned 
position-reference sources J0047-(-5657 (0°.84 separation; 
at ^2000 =00'^47™00.''428805, (52ooo = -h56°57'42"39479 by 
Beasley et al. 2002) and J0042-f5708 (1°.50 separation; 
at a2ooo =00M2"U9."451727, 62000 = +57°08'36."58602 by 
Beasley et al. 2002), switched typically every 10 min- 
utes. A bright calibrator source was also observed ev- 
ery hour at each epoch: J0319-h4130 (= 3C 84) for the 
first epoch and J2232-f 1143 (= CTA 102) for the last five 
epochs. Left-hand circularly polarized signals were dig- 
itized at 2-bit sampling and recorded at a data rate of 
1024 Mbps. Among the total bandwidth of 256 MHz, one 
16-MHz IF channel was assigned to the H2O maser source 
in NGC 281 West and the other 15 IF channels of each 
16-MHz bandwidth were assigned to one of the position- 
reference sources with the VERA digital filter unit (Iguchi 
et al. 2005). The data correlation was performed with the 
Mitaka FX correlator (Chikada et al. 1991). In order to 
achieve a higher spectral resolution for the H2O maser 
lines, only the central 8 MHz of the total 16-MHz band- 
width of the IF channel for the maser lines was split into 
512 spectral points, yielding frequency and velocity reso- 
lutions of 15.625 kHz and 0.21 km s~^, respectively. The 
observed frequencies of the maser lines were converted to 
the radial velocities with respect to the local standard of 
rest (LSR), Vlsr, using a rest frequency of 22.235080 GHz 
for the H2O 6i6-523 transition. The system noise tempera- 
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tures at the zenith were typically 150—300 K for all epochs 
except for the second and fourth epochs, where the system 
temperatures were 400—500 K. The aperture efficiencies 
of the antennas ranged from 45% to 52%. 

Calibration and imaging were performed in a standard 
manner with the NRAO Astronomical Image Processing 
System (AIPS) package. For single-beam data analyses 
to discuss the internal motion of the H2O maser source in 
NGC 281 West, the visibilities of all velocity channels at 
each epoch were phase-referenced to the brightest refer- 
ence maser feature at Vlsr = —31.8 km (feature 4 in 
table 1). After those calibration processes, spectral- line 
image cubes were made using the AIPS task IMAGR. In 
order to search for maser features that were far from the 
reference feature, we first made wide-field maps of 1024 
pixels X 1024 pixels with a pixel size of 1 mas to obtain 
a field of view of ~1", locating the map centers at ex- 
pected positions of the maser components CI, C2 and C3 
(Tofani ct al. 1995). Around the detected maser features, 
narrower-field maps were then made with 512 pixels x 512 
pixels of size 0.05 mas, yielding a field of view of 25.6 mas 
X 25.6 mas, and the positions of these features were de- 
termined with respect to the reference maser feature. The 
synthesized beam has an FWHM beam size of 1.3 mas x 
0.8 mas with a position angle of —43°. RMS noise levels 
per channel were 50—150 mJy. 

Based on our criteria for detection, all maser features 
presented here were detected with a signal-to-noise ra- 
tio higher than 7 over two or more velocity channels. In 
order to avoid the effect of maser structure on position 
measurements, each maser spot was averaged over the de- 
tected velocity channels with the AIPS task SUMIM to 
yield one maser feature, and the feature positions were 
then determined with respect to the reference feature by 
the peaks of Gaussian-model fitting with the AIPS task 
JMFIT. For single-beam multi-epoch data analyses, we 
identified maser features at different epochs as being the 
'same' on the basis that their positions should agree with 
those expected from the proper motions estimated at three 
or more epochs with an accuracy of 1 mas and that their 
radial velocities must be within twice the velocity resolu- 
tion, 0.42 km s~^, from the mean value for all six epochs. 

For dual-beam phase-referencing data analyses to dis- 
cuss the absolute positions and absolute proper mo- 
tions of the H2O maser source, the reduction was done 
in the same manner as the single-beam data reduction 
described above, but using the extragalactic position- 
reference quasars, J0047-I-5657 and J0042-I-5708, as phase- 
reference calibrators instead of the brightest maser feature 
in the single-beam case. In order to correct for instrumen- 
tal phase errors, the dual-beam phase calibration data, 
taken with the artificial noise source in each beam, were 
also included in the dual-beam data reduction. The maps 
obtained from the dual-beam data had low dynamic range 
compared to those of the single-beam data, mainly due to 
residuals of troposphcric zenith delay, and thus only three 
relatively bright maser features (features 4. 9 and 10 in 
table 1) were detected at three or more epochs in the 
dual-beam data reduction. 
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Fig. 1. The observed spectra of the H2O maser source in 
NGC 281 West, obtained from a vector average of the received 
power at the sixth epoch, integrated over the net four-hour 
observation time: (a) the total power spectrum at Mizusawa 
station; (b) the cross power spectrum in Mizusawa-Ogasawara 
baseline. 

3. Results 

3.1. Feature Distribution and Internal Motions 

First, we show the results from single-beam data anal- 
yses to obtain relative positions and internal motions of 
II2O maser features in NGC 281 West. Figures la and lb 
show examples of total and cross power spectra (autocor- 
relation and cross-correlation), respectively. Each spec- 
trum was obtained from a vector average of the received 
power, integrated over the total observation time at the 
sixth epoch (approximately net 4 hours). The peak inten- 
sity occurs at T^lsr ~ —32 km s~^, which is mainly due to 
the reference maser feature. 

Table 1 summarizes the results of both relative and 
absolute proper-motion measurements obtained by linear 
least-squares fits. As listed in table 1, we detected in 
total 10 maser features with radial velocities Vlsr rang- 
ing from —39 to —22 km s~^ at more than two epochs. 
The central velocity of these features is in good agreement 
with the systemic radial velocity obtained from the "'^^CO 
( J = 1 — 0) line emission, Vlsr = —31 km s^^ (Lee & Jung 
2003). Figure 2a shows the distribution of these maser fea- 
tures, which corresponds well with the relative positions 
of two maser components C3 and CI observed with the 
VLA by Tofani et al. (1995). We did not detect any maser 
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Table 1. Results of proper motion measurements. 



ID 




Coordinate 


X 


y 








fJ.y 


Flux 




[km s-i] 


{J2000.0) 


[mas] 


[mas] 


[mas yr"-'] 


[mas yr"-"-] 


[mas yr"-"-] 


[mas yr"-"-] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(C3) 




















1 


-34.1 




222.50 


-55.59 






1.76 (0.03) 


-0.29 (0.04) 


3.9 


2 


-32.5 




-3.39 


1.37 






0.72 (0.14) 


0.09 (0.16) 


1.2 


3 


-32.3 




-5.22 


2.35 






0.30 (0.10) 


-0.07 (0.21) 


1.2 


4 


-31.8 


00'^52™24.=70086 


0.00 


0.00 


-2.87 (0.26)*^ 


-2.78 (0.37)*^ 


0.00 


0.00 


39.8 






+56°33'50."5270 






-2.92 (0.23)*2 


-2.54 (0.35)*2 








5 


-29.8 




137.32 


-138.43 






-1.16 (0.21) 


0.84 (0.09) 


3.5 


6 


-29.6 




101.64 


-89.17 






0.47 (0.04) 


-0.41 (0.06) 


0.6 


7 


-24.3 




83.62 


-14.14 






0.29 (0.07) 


1.48 (0.20) 


0.8 


8 


-22.3 




92.21 


-20.94 






1.33 (0.16) 


1.72 (0.20) 


3.9 


(CI) 




















9 


-39.3 


00'^52m24.=17744 






-4.88 (0.11)*i 


-1.71 (0.46)*i 






5.8 






+56°33'43."2553 






-4.76 (0.28)*2 


-1.44 (0.41)*2 








10 


-32.0 


00'^52™24.n9986 






-2.03 (0.23)*^ 


-4.06 (0.79)*i 






1.1 






+56°33'43."3907 






-1.91 (0.37)*2 


-4.12 (0.69)*2 








4w 










-2.89 (0.18) 


-2.65 (0.26) 








C3m 










-2.43 (0.38) 


-2.23 (0.40) 


0.46 (0.33) 


0.42 (0.31) 




9w 










-4.86 (0.10) 


-1.56 (0.31) 








lOw 










-1.99 (0.20) 


-4.09 (0.52) 








Clm 










-3.43 


-2.83 









(1) The maser feature ID number. Sec text for details. (2) The radial LSR velocity in km averaged for detected epochs. (3) The 
coordinate RA(J2000.0) and DEC(J2000.0) of the absolute position at the first epoch. (4) (5) The best-fit angular offsets in mas at the 
first epoch from the reference maser feature 4, toward the east (RA) and north (DEC) (extrapolated for those not detected at the first 
epoch). (6)(7) The best-fit absolute proper motions in mas yr~^ (with the associated uncertainties in parentheses), obtained using the 
position-reference sources *1: J0047-I-5657 and *2: J0042+5708. (8)(9) The best-fit relative proper motions in the component C3 with 
respect to the reference maser feature 4 (with the associated uncertainties in parentheses). (10) Peak flux in Jy beam"'^ at the sixth epoch. 




Fig. 2. The maps of the maser feature distribution and internal motions: (a) the maser feature distribution, correspond- 
ing well with the relative positions of two maser components C3 and CI observed with the VLA by Tofani et al. 
(1995); (b) the internal motions of features in C3 relative to feature 4, denoted by X at the map origin (see text). 
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Fig. 3. Absolute proper-motion measurements of feature 4; using J0047-I-5657 (open circles) and J0042-I-5708 (open squares) as 
position-reference sources. Note that additional offsets are given for clarity: -1-1.0 mas to the squares in (a) and -1-1.5 mas to the 
circles in (b). The error bars are plotted for the standard deviation a from each linear least-squares fit (sec text). 



emission at the expected position of the component C2. 
With the high spatial resolution of VERA, ^1.2 mas at 
22 GHz, our VLBI map separates the H2O maser emission 
within component C3 into 8 maser features and within CI 
into 2 features. 

In single-beam data analyses, we measured the accu- 
rate relative positions and relative motions of the maser 
features within C3 only, with respect to the brightest ref- 
erence feature 4 in C3 (Vlsr = —31.8 km s~^). Note that 
the other two features located within CI were also used 
for astrometric measurements in the dual-beam phase- 
referencing data analyses, which will be discussed in the 
next section. 

The positions of maser features 1—3 and 5—8 in C3 were 
determined at each epoch relative to the reference feature 
4, and their relative proper motions were also determined 
at detected epochs by linear least-squares fitting. The 
terms x and y in table 1 are the best-fit angular offsets 
in mas at the first epoch from the reference feature 4, 
toward the east (RA) and north (DEC), respectively (ex- 
trapolated from the fits for features 5 and 8, which were 
not detected at the first epoch). In table 1, \ix and [ly 
show the measured relative proper motions with respect 
to the reference feature 4. The associated uncertainties 
in \ix and indicated in parentheses in table 1, were 

estimated from RMS residuals in x and y of 0.1—0.8 mas 
from the linear least-squares fits. 

Figure 2b gives a map of the internal motions of the 
8 maser features in C3, obtained from the best-fit offsets 
(cc, y) at the first epoch with vectors (/.t^, /Xy) of the rel- 
ative proper motions. The reference feature is located at 
the map origin and denoted by X. The best-fit relative 



positions of the features are denoted by filled circles with 
the feature ID numbers in table 1. Note that the magni- 
tude of the vector (/.i^;, /ij,) is magnified to be the value for 
60 yr. A proper motion of 0.5 mas yr~^ is indicated for 
reference at the bottom right corner of figure 2b. This mo- 
tion, 0.5 mas yr~^, corresponds to 6.9 km s~^ at an esti- 
mated distance to NGC 281 of 2.9 kpc (Guetter & Turner 
1997). This map indicates expanding motions of the H2O 
maser features in C3 with velocities from 10 to 20 km . 
As given as C3m in table 1, the unweighted mean of the 
relative proper motions of all 8 maser features in C3 (in- 
cluding feature 4) is (/4, Aty)c3 =(0.46±0.33, 0.42±0.31) 
mas yr~^. The magnitude of the mean relative proper- 
motion vector is 0.62 mas yr~^, which corresponds to 
8.5 km s~^ at a distance of 2.9 kpc. The uncertainties 
associated with the mean were estimated from the stan- 
dard deviation of the relative proper motions, divided by 
\/n—\, where n = 8 is the number of the maser features 
in C3. The standard deviation of these relative proper 
motions is (0.88 ,0.82) mas yr~^, which corresponds to a 
deviation of ~10 km s~^ from the mean motion. 

3.2. Absolute Positions and Absolute Proper Motions 

We now show the results from dual-beam phase- 
referencing data analyses to obtain absolute positions and 
proper motions of the H2O maser features in NGC 281 
West. As shown in table 1, we measured the absolute 
positions and absolute proper motions fix and fiy to- 
ward the east (RA) and north (DEC), respectively, for 
three bright maser features. These measurements were 
performed with respect to two different position-reference 
sources J0047-h5657 and J0042-h5708, independently. For 
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Fig. 4. A map of the absolute proper motions of the maser 
features, indicating the measured proper motions of fea- 
tures 4, 9 and 10 (solid arrows) and the mean proper mo- 
tions of the maser features in CI and C3 (dotted arrows). 
The map origin is at RA(J2000.0)=00'^52'"24.=70086 and 
DEC(J2000.0)= +56°33'50f'5270 (see text). 

the absolute position of each feature, the unweighted 
mean of the two values obtained from these two position- 
reference sources at the first epoch is shown in table 1 
(column 3). The positional difference between the two 
reference cases was ~0.5 mas, likely due to the uncer- 
tainties in the absolute positions of the reference sources: 
0.64 mas for J0047-h5657 and 0.82 mas for J00424-5708 
(Beasley et al. 2002). These systematic errors do not affect 
the precision of the proper motion measurements because 
they only shift the maser feature positions with constant 
offsets. 

Figure 3 shows the absolute proper motion of feature 4. 
The open circles and squares indicate values obtained 
using J0047+5657 and J0042-I-5708 as position-reference 
sources, respectively. Note, however, that additional off- 
sets are given in figure 3 for clarity: -1-1.0 mas to the 
squares (from J00424-5708) in figure 3a; and +1.5 mas 
to the circles (from J0047-(-5657) in figure 3b. The ob- 
servation period of 189 days was not long enough to de- 
tect the annual parallax significantly. Hereafter we as- 
sume the distance to NGC 281 to be 2.9 kpc (Guetter 
& Turner 1997) unless otherwise stated in this paper. In 
figures 3a and 3b, the annual parallax at a distance of 
2.9 kpc was subtracted from the measured absolute posi- 
tion at each epoch. The coordinate of the map origin of 
figures 3 is set to be RA(J2000.0)=00h52'"24."70081 and 
DEC(J2000.0)= -h56°33'50"5274. In figure 3, the error 
bars are plotted for the standard deviation a from each 



linear least-squares fit. Thermal errors in individual po- 
sition measurements with the AIPS task JMFIT ranged 
from 0.010 mas to 0.017 mas (well within the size of each 
circle or square in figure 3), indicating that thermal er- 
rors do not dominate or explain enough the deviation of 
each measurement from the linear fit. We consider that 
the main cause of the deviation from the fit is residuals 
of tropospheric zenith delay, which are difficult to mea- 
sure quantitatively. Wc thus estimate the errors from the 
standard deviations from the fits. 

The absolute proper motion (/ix, ^^Y) of fea- 
ture 4 was measured to be (-2.87±0.26, -2.78±0.37) 
mas yr~^ using J0047-f5657 (denoted by *1) and 
(-2.92zb0.23, -2.54±0.35) mas yr"! using J0042-I-5708 
(denoted by *2). The associated uncertainties in 
and /iy, indicated in parentheses in table 1, were esti- 
mated from RMS residuals from the linear least-squares 
fits to the X and Y offsets, respectively. These RMS 
residuals in X and Y were 0.10 mas and 0.15 mas, re- 
spectively. The error-weighted mean of these two values 
{*1 and *2) is also given in table 1, denoted by 4w, as 
(mx, My) = (-2.89±0.18, -2.65±0.26) mas yr'^. 

The absolute proper motions and /iy of the 

two features in CI (features 9 and 10) were also 
measured in exactly the same manner and are given 
in table 1. The error- weighted mean proper motion 
of each of these features in CI is, as in table 1, 
(mx, My) = (-4.86±0.10, -1.56±0.31) mas yr"! for fea- 
ture 9 and {^JLx, Aiy) = (-l-99±0.20, -4.09±0.52) mas yr'^ 
for feature 10 (denoted by 9w and lOw in table 1, respec- 
tively) . 

Figure 4 gives a vector map of the absolute proper mo- 
tions 4w, 9w and lOw (shown as solid arrows) measured 
for features 4, 9 and 10, respectively. The map origin 
is set to be the absolute position of feature 4 obtained 
at the first epoch: RA(J2000.0)=00''52'"24."70086 and 
DEC(J2000.0)= -b56°33'50"5270. Note that the magni- 
tudes of all vectors {^ix^ fJ-v) are magnified to be the value 
for 500 yr. A proper motion of 2 mas yr~^, which corre- 
sponds to 27 km s~^ at a distance of 2.9 kpc, is indicated 
for reference at the top right corner of figure 4. 

In order to obtain the systemic absolute proper mo- 
tion in C3, we add to the absolute proper motion 
of feature 4 the unweighted-mean vector of the rel- 
ative proper motions of the spots in C3 with re- 
spect to feature 4, which was obtained in the pre- 
vious section to be (/ii;, /Iy)=(0.46±0.33, 0.42±0.31) 
mas yr~^. This addition gives, as denoted by C3ni 
in table 1, the systemic absolute proper motion in C3 
as iiTx, AiV)c3 =(-2.43±0.38, -2.23±0.40) mas yr-^ 
The associated uncertainties were derived from root-sum- 
square calculations of the errors in the mean proper mo- 
tion of C3 (relative to feature 4) and in the absolute proper 
motion of feature 4. This mean motion is plotted in fig- 
ure 4 as a dotted arrow from the mean position of the 
features within C3. 

Since the two features in CI are measured to be moving 
away from each other, as seen in figure 4, we assume an 
expanding internal motion within CI, which is the most 
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Fig. 5. A plot of the mean motion of H2O maser features in C3 and CI relative to the simulated Galactic rotation, superimposed 
on gray-scale and contour maps of the Hi data from Hartmann and Burton (1997). The Hi maps were obtained by integrating over 
Vlsr = -65 to -25 km s-i. The contour peak flux is 75.6 K and the contour levels are 12, 16, 20, 24, 28, 32, 36 K. 



common case of jets or outflows from young stellar objects 
traced by H2O maser motions (e.g., Torrelles et al. 2005). 
We obtain the unweighted mean of the absolute motions 
of these features, as denoted by Clm in table 1, to be 
{fix, AtV)ci =(—3.43, —2.83) mas yr~^. This mean motion 
is plotted in figure 4 as a dotted arrow from the mean 
position of features 9 and 10 in CI. 

If we adopt an estimated distance of 2.2 kpc (Georgelin 
& Georgelin 1976; Roger & Pedlar 1981) or 3.5 kpc 
(Henning et al. 1994) instead of 2.9 kpc in the parallax 
subtraction, the resultant mean proper motions of C3 and 
CI vary only by 10-15%. 

As seen in figure 4, we find that the maser features 
in the two components C3 and CI are systematically 
moving to the southwest. This direction is in agree- 
ment with that of the proper motion of the multiple 
star system HD 5005 (the ionizing source of the Hn re- 
gion NGC 281), which was measured with Hipparcos to 
be [fix, Mr)HD5005 =(-2.95±1.37, -3.22±1.00) mas yr"! 
(Ferryman et al. 1997). The mean proper motion of maser 
features in C3 is (^ix, /iF)c3 =(-2.43±0.38, -2.23±0.40) 
mas yr~^, which agrees with the proper motion of 
HD 5005 within the margin of error of Hipparcos mea- 
surements. While the proper motion of HD 5005 mea- 
sured with Hipparcos over a few years is associated with 
a large uncertainty of 30—50%, we successfully obtained 
the absolute proper motions of the H2O maser features in 
NGC 281 West with a much better precision over a period 
of half a year. 

4. Discussion 

J^.l. Vertical Motion Away from the Galactic Plane 

With the results of our proper-motion measurements, 
we now consider the motion of the H2O maser source in 
NGC 281 West with respect to the Galactic plane. At the 
position of NGC 281 West, the Galactic plane lies almost 
parallel to the east-west (RA) direction with a position 



angle of 90. °2, so that we can regard the motions of the 
maser source in the directions of RA and in DEC as those 
parallel and perpendicular to the Galactic plane, respec- 
tively. 

In this section, wc do not discuss in detail the RA mo- 
tions parallel to the Galactic plane because of their rela- 
tively large uncertainties due to the uncertainties in the 
Galactic rotation. We therefore focus on the motion of the 
maser source in NGC 281 perpendicular to the Galactic 
plane, called the vertical motion hereafter, represented by 
the DEC motion. 

If we assume the distance from the Sun to the Galactic 
center, i?o, to be 8.0 kpc (Rcid 1993), the LSR velocity 
80 of the Galactic rotation at Rq to be 200 km s~^ (e.g., 
Kalirai et al. 2004; Avedisova 2005), and a flat rotation 
curve of the Galaxy, then the expected apparent motion 
of the NGC 281 region due to the Galactic rotation seen 
at the LSR is {hgx, MGy)=(— 2-75, —0.25) mas yr~^ in 
RA and DEC, including the apparent motion in DEC (in 
h) arising from the nonzero galactic latitude of NGC 281 
(6 — — 6.°31). Here the expected radial velocity of the 
NGC 281 at its position (calculated from the flat rotation 
model) is Vi^SR = —31.8 km s~^, which agrees well with 
the observed radial velocities of the NGC 281 West cloud 
and the H2O maser features in the cloud. If we adopt 
the solar motion relative to the LSR based on Hipparcos 
data, corresponding to a solar motion toward the coor- 
dinates a2ooo =251. °51, ^2000 =10.°13 with the velocity 
V =13.34 km s'^ (Dehnen & Binney 1998), we also ob- 
tain the solar-motion effect to the observed motion of the 
NGC 281 region to be (/.tsx, msy) = (0.81, -0.51) mas yr"! 
in RA and DEC. Subtracting these apparent motions due 
to the Galactic rotation and the solar motion from the 
measured absolute proper motion of the H2O maser source 
in NGC 281 West, wc therefore obtain the motion of the 
maser source with respect to the rotation of the Galactic 
Plane. The mean motions of C3 and CI with respect to 
the Galactic rotation arc obtained from the results in the 
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previous section to be {fix, /^y)c3-gr, =(—0.49, —1.47) 
mas yr"i and (fix, a^y)ci-gr =(-1-49, -2.07) mas yr"\ 
respectively. The mean of these two motions is then 
ifJ-x, /^F)sys-GR =(-0.99, -1.77) mas yr~\ which is con- 
sidered to trace the systemic motion with respect to the 
Galactic rotation. Figure 5 depicts this mean motion with 
respect to the simulated Galactic rotation, superimposed 
on grayscale and contour maps based on the H I data from 
Hartmann and Burton (1997). At a distance of 2.9 kpc, a 
proper motion of 1 mas yr^^ corresponds to a transverse 
velocity of 13.7 km s~^, so that the mean vertical motion 
of the maser features is —24.2 km s~^, where the minus 
sign indicates the direction of the vertical motion away 
from the Galactic plane (because of the negative galactic 
latitude b = -6°.31 of NGC 281). 

If we adopt an estimated distance of 2.2 kpc in- 
stead of 2.9 kpc, the resultant mean vertical motion is 
— 18.4 km s~^. If a distance of 3.5 kpc is employed, the 
mean vertical motion is then calculated similarly to be 
-29.4 km s-\ 

If we adopt the lAU standard values Rq —8.5 kpc and 
Go =220 km s"^ (Kerr & Lynden-BcU 1986) instead of 
Rq =8.0 kpc and 8o =200 km s^^ and the solar motion 
relative to the LSR of the velocity V =19.5 km s^^ toward 
"2000 =271. °0, ^2000 =29. °0, then the values vary by up to 
^^10% in the vertical (DEC) motion. 

Our new results from the proper-motion measurements 
clearly demonstrate that the II2O maser features in 
NGG 281 West are systematically moving further away 
from the Galactic plane with a transverse vertical veloc- 
ity of 20—30 km s^^ at its estimated distance of 2.2—3.5 
kpc. 

4-2. Origin of the Motion Away from the Galactic Plane 

In this section, we consider the following three possible 
origins of the detected absolute proper motions of the II2O 
maser features in NGC 281 West, which show a systematic 
motion away from the Galactic plane as described in the 
previous section. 

First, we look at the possibility that the motions arc 
due to outflows of young stellar objects (YSOs) inducing 
the maser emission. Other II2O maser sources show out- 
flow motions of this magnitude (e.g., 10—30 km s^^ in 
W75 reported by Torrelles et al. 2003). Also, the range 
of radial velocities Vi,sR of the maser features (—39 to 
—22 km s^^) are larger than that of the ambient NGC 281 
West cloud (-34 to -26 km s'^ by Lee & Jung 2003) or 
especially of the NW gas clump with maser emission in 
the cloud (-32 to -28 km by Megeath & Wilson 
1997). It may be therefore suggested that the motions 
form in outflows. We regard the internal maser motions 
within C3 or CI as possible outflows in each region, and 
especially the maser features in C3 showed expanding in- 
ternal motions as we have already seen in §3.1 (figure 2). 
However, the internal motions of these maser features in 
C3 are small (~10 km s^^; ~0.8— 0.9 mas yr~^) com- 
pared to the absolute proper motion of the H2O maser 
source, and the central radial velocity of the maser fea- 
tures is in good agreement with that of the entire CO 



molecular cloud complex of the NGC 281 region (includ- 
ing the NGC 281 West cloud), VisR = -31 km s^^ (Lee 
& Jung 2003). The large separation ^^8" between C3 and 
CI in the sky, which is equivalent to a projected separa- 
tion of ~23000 AU (=0.11 pc) at an estimated distance 
of 2.9 kpc, suggests that these two maser components are 
separately produced by two different YSOs. On this scale, 
we only detected the systematic motion in one direction 
and did not detect bipolar jetlike outflows that are often 
traced by H2O maser sources. From these reasons, the de- 
tected systematic mean motion of these two components 
is considered to trace a larger-scale systemic motion of 
the molecular clouds, rather than motions of outflows. It 
is therefore unlikely that the detected systematic motion 
away from the Galactic plane is due to outflows. 

Second, we take into account the possibility that the 
motions are due to the rocket effect of HD 5005, which 
is ionizing the southwestern cloud of NGC 281 (Megeath 
& Wilson 1997). Elmegreen and Lada (1978) suggested 
that the photoevaporation and expansion of the H II re- 
gion NGC 281 is triggering star formation in the NGC 281 
West molecular cloud. Megeath and Wilson (1997) com- 
pared the external pressure of the photoevaporating gas 
to the internal turbulent pressure of the gas clumps in the 
NGC 281 West cloud, and found that the shock velocity is 
low (~1.5 km s~^). Also, we can roughly estimate the ex- 
pansion velocity of the Hll region from an estimated age 
of - 3 Myr of HD 5005 (Henning et al. 1994) and its loca- 
tion 3'A to the east and 3'.8 to the north of the NGC 281 
West cloud (^4.3 pc separation at a distance of 2.9 kpc). 
Dividing the separation ~5' between NGC 281 West and 
the HII region by the age of HD 5005, 3 Myr, we obtain a 
rough estimate of the expansion velocity of 0.1 mas yr""*^, 
which has a negligible effect on the observed systemic mo- 
tion of the NGC 281 region. Therefore, even though the 
rocket effect would accelerate the cloud away from the 
Galactic plane, the effect is unlikely to contribute to the 
detected motions of the H2O maser features significantly. 

Finally, we consider the possibility that the motions are 
due to the superbubble expansion and blowout. Megeath 
et al. (2002) and Megeath et al. (2003) found a broken 
ring of the NGC 281 molecular cloud complex ~ 300 pc 
above the Galactic plane with a ring diameter of 270 pc, 
expanding at a line-of-sight velocity of 22 km s~^ parallel 
to the Galactic plane, which yields a dynamical time of 
the ring of 6 Myr. From the velocity-integrated mass in 
atomic and molecular gas of 3.5 x 10^ Mq and lO^M©, re- 
spectively, they estimated the total kinetic energy of the 
clouds to be 4.5 x 10^^ ergs, which requires multiple super- 
novae (Megeath et al. 2002; Megeath et al. 2003). They 
suggested that the clouds were created from a fragment- 
ing superbubble shell in a superbubble blowout. They 
had no measurement of the motion perpendicular to the 
Galactic plane. From our measurements of the systematic 
motion of the H2O maser features in NGC 281 West away 
from the Galactic plane, it is suggested that the molec- 
ular clouds arc moving at a velocity of 20 — 30 km s~^ 
perpendicular to the Galactic plane, which is comparable 
to the ring expansion velocity of 22 km s~^ parallel to the 
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Galactic plane, thus yielding a comparable kinetic energy 
(~ 4 — 8 X 10^^ ergs) of the clouds for motions perpendicu- 
lar to the Galactic plane as well as parallel. According to 
the numerical model simulation of fragmenting superbub- 
ble blowouts by Mac Low, McCray and Norman (1989), 
the fragments can move at velocities ~ 50— 100 km s~^. 
The observed velocity of the clouds perpendicular to the 
Galactic plane, 20 — 30 km s~^, is lower than the value 
expected from the model, but still consistent with the 
model, taking into account the significant velocity com- 
ponent parallel to the Galactic plane. 

Using the detected vertical velocity v of the region 
away from the Galactic plane, we can also calculate the 
timescale of the explosion for the gas in the NGC 281 re- 
gion to reach the current galactic latitude, assuming that 
its origin is in the Galactic plane. On the simplest assump- 
tion of a ballistic vertical motion, the dynamical timescale 
t is calculated from t = z/v = b/fiy-GRj where z is the 
distance from the Galactic plane, b is the galactic latitude 
of NGC 281 and hy-gr is the derived vertical motion 
of the maser source in NGC 281 West. For b = -6.°31 
and fj-Y-GR = —1.77 mas yr""'^, we obtain 13 Myr as a 
rough estimate of the dynamical timescale of the expan- 
sion, which is comparable in the order of magnitude to the 
dynamical time of 6 Myr by Megeath et al. (2002). This 
timescale of ~ 10 Myr is consistent with the shorter age of 
HD 5005, the ionizing source of the Hll region NGC 281. 
According to earlier studies, the upper limit to the life- 
time of HD 5005, the driving source of the Hll region 
NGC 281, is not much greater than 2 Myr (Elmegreen & 
Lada 1978), and the age of HD 5005 is estimated to be 
~ 3 Myr (Henning et al. 1994). From the difference be- 
tween the estimated timescale of the superbubblc expan- 
sion and the age of HD 5005, it is inferred, as suggested 
by Megeath et al. (2002), that the first OB stars including 
HD 5005 were formed by supernova-driven gas compres- 
sion after the multiple or sequential supernova explosions 
and that the subsequent and ongoing star formation was 
triggered afterward by those first-generation OB star clus- 
ter as suggested by Elmegreen and Lada (1978). 

In conclusion, our measurements of the systemic motion 
of the H2O maser features in NGC 281 West away from the 
Galactic plane provide the most direct evidence that the 
gas in the NGC 281 region and in the HI loop was blown 
out from the Galactic plane, most likely in a superbubble 
driven by multiple or sequential supernova explosions in 
the Galactic plane. 
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